Structure of densified amorphous germanium dioxide 
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Classical molecular dynamics simulations are used to study the structure of densified germanium 
dioxide (GeO?). It is found that the coordination number of germanium changes with increasing 
density (pressure) while pressure released systems do not exhibit a marked change in local structure 
as compared to the virgin system. The structural modification with pressure appears to be stepwise 
and gradually affects long range (through the reduction of the long-range correlations as seen from 
the shift of the first sharp diffraction peak), intermediate range (by angular reduction) and finally 
short range structure (by tetrahedron distorsion). 
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Structural transitions in minerals are known to take 
place under various geological conditions [1,2]. In 
the Earth's interior, silicate and alumino-silicate melts 
change their local structure causing strong density and 
viscosity modifications in magmas [3] and silicon exhibits 
at high pressure a change of its coordination number 
[4]. These results have been mostly obtained with high 
pressure-temperature experiments [5] but also with com- 
puter simulations [6] reproducing extreme conditions. 
Both have inferred the nature of these structural tran- 
sitions, the structure and the phase portrait of the liquid 
state. However, while a majority of studies have been de- 
voted to the silica and silicate chemistry, little has been 
done to elucidate the corresponding behavior in germa- 
nium dioxide (GeO-z) even though this material is a struc- 
tural analog of silica in many respects: both materials 
exhibit at ordinary conditions a tetrahedral local struc- 
ture, they can also exist in a as well as [3 quartz phases 
[7] and the change of germanium coordination from four 
to six also occurs at high pressure [8]. The structure of 
in situ densified or permanently densified GeOi remains 
however controversial. While the global increase with 
pressure of the germanium- oxygen distance in GeOi has 
been related from x-ray diffraction [8] with the conver- 
sion of tetrahedral Ge(4) into octahedral Ge(6), it seems 
that this structural change is reversible as no Ge(6) is 
found in decompressed samples [9] , a situation which docs 
not occur [8] in the decompressed rutilelike c — GeO^ 
The interest in germanium coordination change has been 
also motivated by the observation of the so-called u ger- 
manate anomaly" which corresponds to a maximum in 
density and refractive index when 15 — 16 mol% Na20 
are added into the basic GeOi network former [10] . These 
binary systems have been investigated by various spec- 
troscopic tools [11]- [12] and it is suggested that the in- 
creasing presence of GeOe octahedra (Ge(6)) within the 
network is responsable of the anomaly. On the other 
hand, micro-Raman [13] applied on the same systems 
suggests that the anomaly is due to a massive conversion 
of 4-membered rings into more close-packed ones such as 



3-membered rings, with no, or at least very few, Ge(6) 
present. 

A preliminary task if one wishes to describe density 
induced structural changes in germanates, is first to un- 
derstand how the basic network former GeOi changes 
with dcnsification. While crystalline phases of GeC>2 have 
been studied from numerical calculations [14,15], we are 
not aware of any published result on simulated liquid 
and amorphous GeOi- It is therefore of striking inter- 
est to see what Molecular Dynamics (MD) can tell about 
a certain number of experimental open questions which 
remain at this stage. How does densification affect the 
local structure in the amorphous material ? Which ther- 
modynamical quantity (temperature, Pressure) controls 
mostly the existence of a 6-fold coordinated germanium 
? Does octahedral germanium exists in pressure released 
(permanently densified) Ge02 ? This Letter attempts to 
address some of these basic issues by providing the first 
MD study of amorphous germania. The results show sev- 
eral main features. They demonstrate the reversible na- 
ture of pressure induced changes and no major difference 
with the virgin structure is found upon pressure release, 
while the local coordination number of the germanium in- 
creases smoothly from 4 to 5.5 with applied pressure, a 
result that would be accessible from in situ experiments 
as already realized for silica [16]. Furthermore, several 
regimes of densification can be clearly identified: a first 
regime [P < 1.8 GPa] during which no global change 
in the local structure is found to occur whereas the in- 
crease in the position of the first sharp diffraction peak 
(FSDP) suggests a global reduction of the longer range 
correlations. In the window 1.8 GPa < P < 2.8GPa, a 
sharp decrease of the intertetrahedral bond angle permits 
to densify the structure with no distorsion of the basic 
GeO^/2 tetrahedron. This leads to the buckling of the 
network connected tetrahedra. Finally, for larger pres- 
sures, distorsion of the tetrahadron sets in. As a con- 
sequence, additional constraints appear for P > 3GPa 
which produce a global stiffening of the network [17,18] 

The system consists of 256 germanium and 512 oxy- 
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gen atoms interacting via a Born-Huggins-Mayer type 
potential which has been fitted in the case of GeOi by 
Ocffncr and Elliot to recover the crystalline phases of 
GeOi and its vibrational spectra [14]. The atoms have 
been first confined in a cubic box of length L = 23.044A 
in order to recover the experimental value of the density 
(p g = 3.66 g.cm~ 3 ). After having thermalized the system 
at 3000 K for 20000 time steps (20 ps) the system has 
been cooled to 300 K with a linear cooling schedule at 
a quench rate of 2.5 x 10 12 K.s~ x . Integration has been 
done using a leap-frog Vcrlct algorithm. Various con- 
figurations (positions and velocities) have been saved at 
different temperatures which have been used as starting 
configurations for production runs of 10 5 steps. 



tern at zero pressure was p = 4.5 g.cm~ 



1.25p g . The 




FIG. 1. Partial structure factors Sij(Q) of vitreous germa- 
nia at 300 K from MD simulations (solid line) at p = p g 
compared to experimental results (circles) from [11]. 

Dcnsification has been realized reproducing experi- 
mental conditions, i.e. starting from an initial config- 
uration at T = 300 K and p g and increasing the density 
during 10 4 time steps. At the density p g , the glass tran- 
sition temperature was about 1660 K, and as expected, 
the T g shifts to the higher temperatures with increasing 
density (1894 K at p = l-lp g ). The sample data of Price 
and co-workers [12] has been almost recovered, i.e. af- 
ter MD pressure release from 15.16GPa determined from 
the virial [19]. The final density of the decompressed sys- 



corresponding configuration (termed in the following as 
11 permanently densified system") was taken for compari- 
son with the experimental study [20] of permanently den- 
sified Ge0 2 

The results for the partial structure factors Sij (Q) are 
displayed in Fig. 1 at p g — 3.66 g.cmT 7 ' which show a fair 
agreement with the experimental results obtained from a 
combination of x-ray data [121. 
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FIG. 2. Calculated partial pair correlation functions in vit- 
reous germania at 300K at ordinary density p g — 3.66 g.cm~ 3 
(solid line) and in a permanently densified system 
(p = 4.5 g.cm~ , dashed line). The insert represent the corre- 
sponding running coordination number nij(R) together with 
results from a system under 16.6 GPa pressure (dots). 

One has however to note that these measurements have 
been refined but without any complete resolution of all 
the partials [20]. The principal features in the struc- 
ture factors of GeOi are the peaks occuring at Qp = 
1.5 — 2 AT 1 (the first sharp diffraction peak, FSDP, cor- 
responding to a correlation length L\ — 2-k/Qp ~ 4.lA). 
Peaks at 2.62" 1 and 4l _1 appearing in the partial struc- 
ture factors of the Ge-0 and 0-0 can be associated 
with chemical short range order (CRSO) and topological 
short-range order as derived from scaling considerations 
[21]. For the former, the similar height in goeGe and goo 
but with opposite sign suggests the CRSO nature of the 
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network. The evolution of the total scattering function 
and the position of the FSDP with respect to pressure is 
discussed below. 

The calculation of the pair correlation functions per- 
mits to extract the simulated bond lengths which also 
agree with experimental findings: doe-Ge — 3.321, 
dce-o — 1-72 A and do-o — 2.811 (to be respectively 
compared to the experimental values [11] of 3.16±0.03l, 
1.73 ± 0.031 and 2.83 ± 0.051). Weak changes are found 
when comparing the virgin and the permanently densified 
systems. The most striking difference is found in the Ge- 
Ge correlator which differs even at low distances whereas 
the loss of peak heights in goo produces a slight shift 
in the corresponding running coordination number noo- 
The number of oxygen neighbours around an germanium 
atom remains however unaffected by pressure release. 
More dramatic is the change in structure in the pressur- 
ized system at 16.6 GPa (dots in the inserts of Fig. 2) 
which show substantial differences in the bond distances 
(doe~Ge = 3.25 1, d Ge _ = 1-75 1 and d -o = 2.56 1) 
and an increase of the number of oxygen neighbours in 
the vicinity of a Ge atom (about 5 in the insert represent- 
ing nceO, see also Fig. 3). Evidence of supplementary 
atoms in the first shell surrounding a central Ge0 4 / 2 unit 
is also provided by the increase of nceGe- 




FIG. 3. Left axis: The distorion parameter of a regular 
GeO/4/2 as a function of applied pressure(open circles). For 
comparison, the same parameter for SiOi (filled circles). The 
dashed horizontal line represents the value -^3/8. The right 
axis represents the Ge-0 coordination number (solid line). 

Using the present simulation, it is now possible to an- 
alyze in more depth the low pressure behaviour. Focus 
is made on the distorsion of the tetrahedron parameter 
defined by 6 = dx-o/do-o [X=Si,Ge] which provides 
a direct measure of the effect of the pressure on the lo- 
cal structure of the network. For an ideal tetrahedron, 
S = y/3/8. The latter quantity is of central interest when 
studying the flexibility of the glass under pressure [18] 
and the rigid unit modes in the context of pressure in- 
duced rigidity [22,23]. In Fig. 3 (left axis) are represented 
the variation of the tetrahedron parameter S with pres- 



sure for both silica [24] and germania. At low pressure, 5 
remains almost constant suggesting that the tetrahcdral 
environment is preserved, slightly higher however than 
the value of a perfect tetrahedron. For P > 2.8 GPa, 
there is increasing distorsion of the tetrahedron in ger- 
mania. The way of distorsion appears however to be 
radically different as compared to the silica system. In 
the former it is found a stepwise increase (most noticable 
from the jump of 5 at around 3 GPa), in contrast with 
the more or less smooth increase of 5 for the latter. In 
the insert of Fig. 3 are shown the variation of the mean 
bond angles O-Ge-0 and Ge-O-Ge with respect to pres- 
sure. For both the virgin and the permanently densified 
system, the O-Ge-0 bond angle peaks at around 109°. 
As pressure is increased, the intratetrahedral bond angle 
decreases since the Ge-0 bond distances increases. The 
constant value of the O-Ge-0 bond angle at low pres- 
sure correlates of course with the absence of distorsion 
of the basic tetrahedron. The main feature provided by 
the angular analysis comes from the variation of the in- 
tertetrahedral mean bond angle Ge-O-Ge with pressure, 
which exhibits a sharp drop at around 1.8 GPa, from 
158° to 135° followed by a stabilization at around 140°. 
Densification first applies on the angles connecting the 
tetrahedra and preserves the latter. If pressure keeps in- 
creasing, the distorsion of the tetrahedra sets in, which 
in turn stabilizes more or less the Ge-O-Ge bond angle. 
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FIG. 4. Scattering factor S(Q) for different applied pres- 
sures: 0,3.0, 5.76, 8.92 and 16.6GPa. The lower curve is the 
virgin system at OGPa, compared to experiments [11]. The 
insert shows the simulated position in Qp of the FSDP with 
pressure (solid line). Experimental data are from [23] (open 
circles) and [9] (filled circles). 

Further interpretation is provided from the variation 
of the position of the FSDP with respect to pressure. In 
Fig. 4 are represented the scattering functions for dif- 
ferent applied pressures which show a global broadening 
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of the peak at 4 A" 1 whereas the position of the FSDP 
at 1.5 A -1 is shifted to the higher values in Q, already 
at very low pressure and even before the angular drop at 
P = 1.8 GPa. With increasing pressure, the calculated 
FSDP broadens and becomes less intense, as currently 
observed from X-ray or neutron diffraction [11]. It is 
worth to note that the peak at 2.5 A -1 obtained in the 
present simulation is only weakly observable in the ex- 
periments from neutron diffraction [9] displayed in Fig. 
4. However, the simulated double peak distribution be- 
tween 1.5 and 2.5 A -1 has been observed by different 
authors [26] (see also discussion in [12]). The evolution 
with pressure of the position Qp of the FSDP is repre- 
sented in the insert of Fig. 4. Both experiments and sim- 
ulation show that Qp already increases at low pressures 
and then stabilizes at around Q — 1.7 A -1 (~ 2 A -1 in 
the simulation). This suggests that intermediate range 
order is immediately affected by the densification and 
then remain unaffected with further densification. 

In conclusion, we have shown that simulated GeOi un- 
der pressure shows several main features with applied 
pressure and pressure release: i) a global increase of the 
number of oxygen neighbours in the vicinity of a germa- 
nium atom, ii) a stepwise change in the local structure 
with applied pressure, made of a reduction of long range 
correlation (seen from the position of the FSDP), a sharp 
reduction of the intertetrahedral bond angle and then a 
progressive distorsion of the Ge0 4 / 2 tetrahedron and iii) 
no noticable change in local structure between a virgin 
and a permanently densified system. This clearly draws 
the following picture: Pressure applies on different length 
scales. With increasing magnitude, densification is real- 
ized by a successive deformation of long range structure, 
intermediate (angular) and finally short range structure 
(tctrahcdral) . 
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